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摘 要 : 为 揭示 锡 林 河 流域 上 游 燕 散 发 CET) 变化 规律 及 其 时 空 分 布 特征 ， 采 用 Landsat 8 数据 、 气 象 观测 
数据 和 SEBAL 模型 估算 2015—2017 年 共 12 期 植物 生长 季 的 蒸 散 量 ， 并 分 析 其 时 空格 局 及 影响 因素 。 结 
果 表 明 : O 研究 区 模型 反 演 参 数值 与 实测 值 拟 合 效果 良好 ， 地 表 净 辐射 《R,)、 土 壤 热 通 量 (G)、ET 的 决 
定 系数 (CR?) 分 别 为 0.80、0.65 和 0.86; © 同一 时 段 的 Ra A ET 空间 分 布 具有 较 高 的 一 致 性 ， 整 体 呈 现 东 
高 西 低 的 趋势 ， 具 有 显著 的 时 空 分 异 特征 。4 月 日 ET 值 在 0-4.67 mm, 5—7 AA ET 值 整体 呈 增 加 趋势 ， 
最 大 可 达 10.32 mm, 9—10 AH ET 值 减少 至 0-3.34 mm; © 不 同 土地 利用 类 型 蒸 散 量 依次 为 : 水 体 > 沼 泽 
地 > 耕地 > 草地 > 沙 地 > 农村 居民 点 ，ET 值 与 日 平均 气温 《7T)、 饱 和 水 气压 差 (VPD)、R 呈正 相关 ， 与 空气 
WERE CRA) 呈 负 相关 ， 对 Ri, 的 变化 最 为 敏感 。 

关键 词 : SEBAL 模型 ， 莹 散发 ， 时空 特征 ; 锡 林 河 流域 


流域 燕 散 发 〈evapotranspiration，Z7) 包括 土壤 蒸发 、 水 面 蒸 发 和 植被 蔡 腾 ， 是 地 表 水 
热平衡 的 重要 参量 品 ， 可 为 植物 的 生长 发 育 提供 必要 的 水 分 和 能 量 ， 对 干旱 半 干 旱地 区 内 陆 
河流 域 生态 水 文 循环 过 程 起 着 至 关 重要 的 作用 外。Rosenberg 等 中 指出 ， 地 表 降 水 的 70% 以 
蒸 散 发 形式 回 到 大 气 中 ， 在 干旱 区 该 比例 可 达 90%。 因 此 ， 准 确 估算 流域 蒸 散 发 的 时 空 变 
化 ， 对 于 了 解 干 旱 半 干旱 地 区 水 文 循环 过 程 和 区 域 气候 变化 ,实现 生态 、 水 文 、 社 会 绿色 协 
调 发 展 具有 重要 意义 。 

从 1802 年 道 尔 顿 (Dalton) 提 出 燕 发 量 与 水 汽 压 差 成 比例 关系 申 ， 到 20 世纪 50 年 代 估算 
潜在 蒸发 的 Penman 公式 的 提出 马 ， 再 到 1965 年 Monteith 引入 冠 层 阻抗 的 概念 进一步 完善 
了 模型 中 ， 人 们 在 为 实际 燕 散 量 的 理论 估算 进行 了 长 期 不 懈 的 努力 。 但 传统 估算 藻 散 发 的 方 
法 大 多 基于 观测 站 点 (如 波 文 比 系统 等 )， 局 限于 单 点 尺 BE, 很 难 在 区 域 尺度 扩展 , 空间 代表 
ERAN, 随 着 遥感 技术 的 发 展 , 极 大 地 方便 了 大 尺度 非 均 匀 下 垫 面 的 地 表 特 征 参数 的 获取 ， 
可 以 精确 估算 地 表 实 际 燕 散 量 的 一 系列 遥感 模型 应 运 而 生 四 。 例 如 SEBALP'. METRIC”, 
SEBSI , 其中，SEBAL 模型 是 基于 地 表 能 量 平衡 原理 来 反 演 燕 散发 的 方法 ， 模 型 假设 区 
域 存在 冷 热 极 限 点 ， 很 好 的 考虑 了 干旱 地 区 空间 异 质 性 52， 不 仅 在 湿地 、 耕 地 、 居 民 点 等 
土地 利用 类 型 有 较 好 的 适用 性 ， 而 且 对 水 体 藻 散 的 反 演 有 较 高 的 精度 中 时 。 刘 曼 晴 等 号 利 
] SEBAL 模型 反 演 了 水 面 蒸 发 且 以 此 为 估 测 值 ， 与 燕 发 严实 测 值 有 好 的 拟 合 效果 。 连 晋 嫌 
等 09 基 于 SEBAL-METRIC 模型 估算 了 夏季 黑河 中 游 样 带 尺 度 不 同 土地 覆盖 类 型 蒸 散 量 , 结 
果 表 明 : 在 荒漠 区 与 绿洲 区 分 别 选 取 “ 热 点 "可 有 效 提高 模型 估算 精度 。Lee 等 "利用 改进 的 
SEBAL 对 韩国 的 蒸 散 量 进行 了 计算 ， 得 出 平原 地 区 的 蒸 散 量 高 于 高 原 地 区 的 结论 。 

锡 林 河 流域 是 我 国 的 内 陆 河 流域 ,由 于 水 汽 补给 不 足 ， 加 之 不 合理 人 类 活动 的 影响 ， 导 
致 极端 水 文 事件 频 发 、 径 流 调 节能 力 降低 、 草 原 植被 矮小 且 多 样 性 降低 、 河 谷 湿 地 萎缩 增多 ， 
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流域 生态 水 文 功能 退化 等 现象 日 趋 严 重 09。 因 此 亟 需 开 展 内 陆 河流 域 蒸 散 时 空 变化 特征 研 
究 ， 掌 握 流 域 生 态 恢复 的 耗 水 规律 ， 以 期 为 该 区 域 植 被 对 水 热 分 配 响应 提供 科学 依据 。 
1 研究 区 概况 及 数据 
1.1 研究 区 概况 

锡 林 河 流域 位 于 内 蒙古 中 部 锡林郭勒 盟 ,发源 于 赤峰 市 克 什 克 腾 旗 宝 尔 图 山 , 在 贝 力克 
牧场 转向 西北 流 经 锡林浩特 市 ,最终 注入 查 干 诺 尔 沼 泽 地 自然 消失 。 流 域内 由 小 型 湿地 、 滩 
地 、 沼泽 地 、 河渠 等 水 域 组 成 , 是 典型 的 干旱 半 干 旱 区 内 陆 河 流域 。 好 来 吐 郭 勒 、 好 来 郭 勒 、 
呼 斯 特 河 是 其 最 主要 的 3 条 支流 09， 主 河 全 长 175 km (图 1)， 锡 林 河 河源 至 库 尼 苏 曼 为 流 
域 上 游 (43926'-44*08'N. 116°02’~117°12'E), LEJXBIX] 12 4) 828 Cf] SRA ETS (Stipa 
grandis). Æ% (Leymus chinensis), KEZ (Bromus inermis). 等 优势 物种 。 土 壤 类 型 主 
要 以 黑 钙 士 、 栗 钙 圭 为 主 ， 土 壤 有 机 质 含量 较 高 ， 腐 殖 质 层 平均 37cml”"。 该 区 域 气 候 属 大 
陆 性 温带 半 干 旱 气 候 ， 流 域 多 年 平均 降水 量 约 278.9mm， 其 中 6—9 月 降水 量 占 年 总 降水 量 
的 80% 以 上 ， 年 平均 气温 为 2.8CPC。 
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1 研究 区 示意 图 2 土地 利用 类 型 图 
Fig.1 Location of study area Fig.2 Map of land use types 


1.3 遥感 数据 

土地 利用 类 型 数据 提取 于 2015—2017 年 的 Landsat 8 影像 , 根据 我 国 《 土 地 利用 现状 分 
类 标准 》 将 研究 区 共 分 为 7 类 (图 2)， 包 括 水 体 (河流 和 湖泊 )、 沙 地 、 耕 地 、 沼 泽 地 、 高 
履 盖 草地 、 低 履 盖 草地 、 农 村 居民 点 。 

本 文 所 采用 的 遥感 数据 为 Landsat8 OLI 数据 。 在 数据 选择 方面 ， 由 于 生长 季 植 被 覆盖 
度 高 ， 基 于 SEBAL 模型 反 演 的 蒸 散 量 更 加 准确 ， 根 据 研究 区 已 有 卫星 数据 的 影像 质量 、 成 
像 时 间 等 信息 , 选取 锡 林 河 流域 上 游 2015—2017 年 生长 季 的 12 W (2015 4E 4.5. 7. 9 月 ， 
2016 4E 5. 6. 9H, 2017 4E 4. 5. 6. 7. 10 HO 无 云 或 少 云 ， 清 晰 度 高 的 遥感 影像 ， 影 像 
数据 取 自 美国 地 质 调查 局 (USGS) 网 站 (http://earthexplorer.usgs.gov/) 以 及 地 理 空 间 数 据 云 网 
站 (http://www.gscloud.cn/)。 

1.4 气象 数据 及 其 他 辅助 数据 

选取 研究 区 2015 一 2017 年 各 年 4 一 10 月 石门 气象 站 的 气象 数据 ,包括 气温 、 相 对 湿度 、 
风速 、 日 照 时 数 、 降 水 和 土壤 温度 等 。 数 据 采 集 频 次 为 10 Hz， 每 30 min 记录 1 次 ， 通 过 数 
据 采 集 器 CR1000(CR1000, Campbell Scientific) 储 存 ， 仪 器 型 号 和 安装 高 度 见 表 1。 

为 提高 模型 的 估算 精度 ， 在 SEBAL 计算 过 程 中 输入 研究 区 的 气温 和 风速 数据 。 本 文 利 
用 研究 区 附近 的 4 个 国家 站 2m 处 气温 (7T,) 和 风速 (W), 通过 克 里 金 插值 法 (Kriging method) 
得 到 与 遥感 数据 空间 分 辩 率 相同 的 栅 格 数据 。 

使 用 的 DEM 数据 为 地 理 空 间 数 据 云 网 站 (http:/www.gscloud.cn/) 下载 的 GDEMV2 30 m 
分 辩 率 数字 高 程 数 据 。 


表 1 气象 站 传感器 基本 信息 


Table 1 Basic information of micrometeorological sensors 


产品 名 称 型 号 安装 高 度 /m 生产 商 
空气 温 湿度 HMP155A 2 Campbell Scientific Inc., USA 
四 分 量 净 辐 射 仪 NR-LITE 2 Campbell Scientific Inc., USA 
风 杯 式 风速 仪 034B 2, 3.5 Met One Instruments Inc., USA 
土壤 热 通 量 板 HFP01 -0.1 Husker Flux Inc., USA 
雨量 简 TEM525MM 2 Texas Electronics., USA 
土壤 三 参数 仪 CS655 -0.1, -0.2, -0.4, -0.8 Campbell Scientific Inc., USA 


2 研究 方法 
2.1 SEBAL 蒸 散 模型 及 其 参数 计算 
SEBAL 模型 是 Bastiaanssen 等 "提出 的 单 层 计 算 模型 ， 利 用 能 量 平衡 方程 计算 燕 发 瞬 
时 潜 热 通 量 。 然 后 对 图 像 的 每 个 像素 点 计算 瞬时 燕 发 和 日 蒸发 。 
po 6-—H (1) 
sh: 4 是 水 的 汽化 潜 热 (W.m”*.mm!); ET ERRE (mm); R, 为 净 辐 射 量 (W.m”); G 
为 土壤 热 通 量 (W.m”); H JERBE COW-m?-mm D. 
(OD 净 辐 射 通 量 Rd 计算 公式 : 
Rn = (1— oa)Rs, + Ruy — Rir — (1 — €0) Ru (2) 
式 中 : a 为 地 表 反 照 率 ; 人 t UK EH ESAE (Wm); e 
太阳 入 射 的 长 波 辐射 (W.m”); RL 为 地 表 发 射 的 长 波 辐射 (W.m”); e 为 地 表 比 辐射 率 ， 介 于 
0-1 之 间 。 
(2) 地 表 反 照 率 计算 公式 为 : 
g= ee Sean a (3) 


Tsw 


RP: atoa 为 宽带 行星 大 气 顶 反射 率 ; wapath_radiance 是 大 气 路 径 反 照 率 , (A 0.025 — 0.0471; 
Tow 是 大 气 单 向 透射 率 ， equ 


490754 2x10? xZ (4) 
(3) 入 射 的 短波 辐射 TN us 公式 : 
Ry = G4, X Cos0 X d, X Tsw (5) 


式 中 : Gwe 是 太阳 常数 , 取 1367 (Wm?) 9 为 太阳 天 顶 角 (从 卫星 图 像 的 MTL 文件 中 获取 ); 
d, JI BHEN. 
(4) 入 射 的 长 波 辐射 RL》 计算 公式 : 
Ru = &oT; (6) 
RP: e 是 大 气 比 辐射 率 ，G 为 斯 蒂 芬 - 玻 尔 兹 曼 常数 ， 取 5.67x103 W-m?- K^; 7 为 空气 温 
度 。 


£a = 1.08( 一 Inrsw)0265 (7) 
(5) 地 表 长 波 辐射 (RLt): 
Rix = sc0074 (8) 
WP: go 为 地 表 比 辐射 率 ; o IER T 为 地 表 温 度 〈 开 )。 
£o = 1.009 + 0.047 - In(NDVI) (9) 
(6) 土壤 热 通 量 (G): 利用 植被 指数 、 地 表 净 辐射 与 其 关系 , oe PEU TOROS i (E P. 
c= Ts (10 


Rn (0.0038a+0.0074a2)(1—0.98NDVI4) 


(pxCpxdT ) 


H = (11) 


Tah 

式 中 :p 为 空气 密度 , 标准 状态 下 取 1.293 kg-m^; C 空气 定 压 比 热 , 通 常 取 1004 (P kg  -K 1); 
dT 是 高 度 为 Z1 和 Z 处 的 温度 之 差 ，r 为 空气 动力 学 阻力 (S-m") 。 

为 求 得 dT 需要 从 影像 中 选择 干 湿 像 元 : “热点 ”是 植被 覆盖 很 少 的 盐碱地 , 温度 较 高 。“ 冷 
点 ”处 水 分 充足 ， 植 被 茂盛 ， 地 表 温 度 很 低 ， 通 过 计算 “热点 ”和 “ 冷 点 ”就 可 以 得 到 所 有 像 元 
H3 aT fü. 
22 蒸发 比 及 日 蒸 散发 量 

由 于 卫星 过 境 得 到 的 是 瞬时 影像 因此 , 需要 转换 得 到 日 通 量 . SEBAL 模型 中 假设 24h 
ZA ARAL EAT FE AY o 


蒸发 比 人 : 
AET 
Ainst= RG —^24 (12) 
AP: A24 是 一 天 24h 内 的 蒸 散 发 比 。 
PRU ETa: 
_ Rp24X^24X86400 
ET 24 = [2.501—0.002361x(T,—273.15)]x106 (13) 
AP: Rut RAR ENE Wma: T, heil E CK) 
2.3 FAO Penman-Monteith 模型 
其 计算 公式 为 pa 
ET, = 04084(Rs -G) 7.77 uz (65-63) (4) 


pty(1+0.34u2) 
ET = K, x ET, (15) 

RP: E70 为 潜在 蒸 散 发 (mm-d"); ET ASA (mm.d 0: 及 为 作物 系数 ，9 为 
饱和 水 汽 压 曲线 斜 率 CkPa-?C 0: R 为 冠 层 表面 净 辐 射 C(MJ.m .d7); 了 为 日 平均 气温 CC ); 
G 为 土壤 热 通 量 (MI-m?-d'); u Jy 2m 高 度 处 风速 Om:s Di e, 为 饱和 水 汽 压 (kPa): e, 
为 实际 水 汽 压 (kPa); y 为 湿度 计 常 数 (kPa C” )。 
2.4 模型 适用 性 评价 

为 评价 SEBAL 模型 的 适用 性 ,采取 相对 均 方差 (RMSE)、 一 致 性 指数 (d)、 平均 偏 差 
(bias)、 和 Nash-Sutcliffe MEAL (E) 4 个 统计 量 进行 评价 ， 计 算 公 式 如 下 29。 


i 

RMSE = G N(R 0)*f (16) 
-1_ | Zo — C 

d B 1 [Be xl S 

bias = +X}, (0; — Pi) (18) 

E = 1 — [Xi.4(0; - Pi)? / (0; - 0)7] (19) 


2.5 敏感 性 分 析 
敏感 性 分 析 是 从 定量 角度 分 析 相 关 和 气象 因 子 发 生变 化 时 ， 对 ET 的 变化 产生 影响 的 定量 
BBO, HG ONEN ET 变化 率 与 气象 因子 变化 率 的 比值 。 


(200 


= Su . ET VQ 
S= limav-o( AQ/Q ) dQ ET 


: 5 为 敏感 系数 ; AE7 为 实际 蒸 散 量 的 变化 量 ; Q 和 AQ 分 别 为 气象 因子 及 其 变化 量 。 


式 
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3 结果 与 分 析 
3.1 SEBAL 模型 反 演 结果 验证 

为 验证 SEBAL 模型 反 演 精度 ， 本 文 基于 石门 湿地 气象 站 气象 数据 (石门 气象 站 的 建站 
时 间 为 2015 年 11 月 ， 所 以 对 模型 的 验证 从 2016 年 生长 季 开 始 )， 以 及 FAO 56 P-M 计算 的 
蒸 散 量 与 卫星 过 境 时 刻 SEBAL RIR ET, Ra G 值 进行 验证 对 比 。 作 物 系数 K. EE 
长 初期 、 生 长 中 期 和 生长 后 期 取 值 分 别 为 Kon=0.4、Kemia=0.93 和 Keena=0.8 0 

由 实测 值 与 反 演 值 的 回归 分 析 (图 3) 可 知 : ET 的 模拟 值 与 实测 值 的 分 布 较为 集中 ， 
大 多 分 布 在 1:1 直线 附近 ， 决 定 系 数 RR 最 高 ， 为 0.86，R 和 G 均 在 1:1 线 之 上 ， 反 演 值 大 
FSW, ARH, Ry 决定 系数 RR 为 0.80，G 的 模拟 值 较 分 散 且 波动 相对 较 大 ，R" 最 
小 为 0.65; 由 统计 参数 C 2) 可 知 : RMSE 可 衡量 模拟 值 与 实测 值 之 间 的 偏差 , ET 的 RMSE 
为 0.95， 表 明 模 拟 值 与 实测 值 较 接 近 ， 模 拟 结果 最 好 ; d 值 越 接 近 1， 表 明 模 型 模拟 效果 越 
If, Ray Gy ET AHA 0.99, 0.99 和 0.96， 均 趋 近 于 1， 模 拟 效 果 较 好 ; bias 反映 实测 值 与 
模拟 值 的 偏差 越 接近 0， 精 度 越 高 ， 而 Nash-Sutcliffe 效率 指数 (E) 越 接 近 1， 表 示 模 型 可 
信 度 越 高 。 综 上 所 述 ，SEBAL 模型 对 各 个 参数 的 模拟 精度 较 高 ， 因 此 ，SEBAL 模型 在 锡 林 
河流 域 上 游 具 有 较 好 的 适用 性 。 


y = 1.0669x + 67.035 y=0.8953x + 0.252 m y 2 0.82x + 15.821 
^ R22 0.7984 8 R?= 0.862 R2- 0.6515 
一 70 
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2 
at 1 Od 30 Q 
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3 实测 值 与 模型 反 演 值 对 比 
Fig. 3 Comparison of measured values and model inversion values 
表 2 SEBAL 模型 适用 性 分 析 
Table2 Adaptability analysis of SEBAL model 


参数 RJ(W-m?) G/(W-m?) ET/imm 
相对 均 方差 (RMSE) 47.21 10.77 0.95 
一 致 性 指数 (q) 0.99 0.99 0.96 
平均 偏差 (bias) -67.88 -7.14 0.48 
Nash-Sutcliffe 效率 指数 (E) 0.94 0.95 0.85 


SEBAL 模 型 估算 日 ET 值 与 FAO 56 P-M 计 算 值 对 比 结 果 见 图 4， 由 图 4 可 知 : 2016—2017 
年 石门 湿地 生长 季 SEBAL 模 型 估计 值 与 P-M 模 型 的 计算 值 变化 趋势 基本 一 致 , 二 者 起 伏 大 致 
相同 。 营 散 量 的 季节 变化 较为 明显 ， 研 究 区 2016 一 2017 年 4 一 5 月 蒸 散 量 处 于 较 低 水 平 ， 平 均 
值 为 4.81 mm.d ， 从 6 月 开始 呈 增 加 趋势 ，6 一 8 月 苔 散 量 平均 值 为 5.76 mm-d', 9H RAEN 

台 减 弱 ，9% 一 10 月 蒸 散 量 平均 值 为 2.89 mm.d  。 在 出 现 降雨 的 时 间 ， 大 部 分 BE7 都 会 对 应 出 现 
个 波 谷 ， 降 雨 事件 结束 之 后 会 出 现 一 个 波峰 。 
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加 降水 ”一 一 P-M 计 算 值 e 遥感 反 演 值 


4 2016—2017 年 石门 湿地 降水 及 SEBAL 与 P-M 的 计算 值 
Fig. 4 Precipitation, SEBAL and P-M calculation value of ShiMen Wetland from 2016 to 2017 
3.2 R 和 ET 的 时 空格 局 
利用 SEBAL 模型 估算 研究 区 2015—017 年 4—10 月 的 日 Rs 值 , Rs, 时空 分 布 格局 如 图 5 
所 示 。 
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5 2015—2017 年 SEBAL 估算 日 R 时 空 分 布 


202007.00037v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


Fig. $ Spatio-temporal distribution of SEBAL estimated daily net radiation fluxes of 2014 to 2016 


由 图 5 可 以 看 出 ， 锡 林 河 流域 上 游 的 地 表 净 辐射 具有 明显 的 季节 性 变化 ， 总 体 R, 值 在 


100-780 W.m”， 最 大 值 和 最 小 值 相差 很 大 ， 空 间 上 呈现 东 高 


西 低 的 状态 。 生 长 季 内 河谷 湿 


地 的 R, 一 直 保 持 较 高 的 水 平 ,， 在 400-750 W.m?， 在 相同 气候 条 件 下 ,不 同 土地 利用 状况 直 
的 分 配 ， 与 其 他 土地 利用 情况 相 比 ， 河 谷 湿地 水 分 充足 ， 植 被 覆盖 度 高 ， 热 
容量 大 ， 而 地 表 温 度 和 地 表 反 照 率 低 ， 其 对 太阳 辐射 的 吸收 能 力 强 ， 反 射 能 力 弱 ， 相 应 的 地 


接 影 响 地 表 能 量 


表 净 辐射 值 较 大 , 所 以 在 生长 季 内 其 净 辐 射 通 量 值 远 高 于 其 他 土地 利 月 


日 类 型 ; 草地 和 沙 地 的 


R, 在 4 月 保持 在 100-300 W-m*, 5—7 ABA SIM, ZF 300-650 W-m?, 说明 此 时 植被 生命 


活动 旺盛 ， 太 阳 辐 射 是 主要 的 能 量 来 源 。9- 一 10 月 R, 值 整体 降低 ， 为 120-520 W-m?. 


研究 区 基于 


的 变化 趋势 ， 在 5 一 7 月 随 着 日 照 时 间 和 地 表 净 辐射 的 增加 ， 营 散 量 也 随 之 增 大 。 
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研究 区 内 不 同 土地 利用 / 履 被 类 型 及 不 


SEBAL 模型 影响 估算 的 2015—2017 年 日 燕 散 空间 分 布 情况 如 图 6 所 示 。 
对 比 图 5 和 图 6 可 以 看 出 ,地 表 净 辐射 越 高 ， 蒸 散 量 也 越 大 ， 地 表 净 辐射 与 蒸 散发 呈现 相同 
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图 6 2015—2017 年 SEBAL 估算 日 ET 时 空 分 布 
Fig. 6 Spatio-temporal distribution of SEBALestimated daily evapotranspiration of 2014 to 2016 


同时 间 的 莹 散发 差异 明显 ， 整 体 上 来 看 
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呈现 东 高 


西 低 的 趋势 ， 东 部 为 锡 林 河 流域 发 源 地 ,河谷 湿地 植被 高 度 和 履 盖 度 都 较 高 ， 以 植被 蒸腾 为 


E, Ak, RAEAN. 
BIO I vuU, HALERE Rt 2 


H 


mm.d ， 主 要 是 湖泊 和 水 库 ; 


部 草原 植被 覆盖 度 低 ， 土 壤 表 层 含水 量 少 ， 地 下 水 对 地 表 蒸 


行 统 计 ， 研 究 区 最 大 ET 可 达 10.32 


河谷 湿地 蒸 散 无 明显 的 地 形 差异 ， 日 蒸 散 量 空间 分 布 较 均 一 ， 
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蒸 散 量 一直 保 持 较 高 的 水 平 ， 在 2.42-10.32 mm-d'; 耕地 在 作物 生长 季 5—9 FRAC KOE 
仅 次 于 湿地 ， 在 1.82~8 mm-d'; 图 7 中 部 分 农村 居民 点 日 蒸 散 量 也 较 高 ， 这 是 由 于 一 些 农 
村 居民 点 附近 种 植 蔬菜 等 作物 , 有 旦 农村 居民 点 一 般 与 农田 毗邻, 存在 较 多 与 农田 的 混合 像 元 ; 
沙丘 裸 地 土壤 含水 量 低 ， 植 被 蒸腾 作用 小 ， 蒸 散 量 少 , 在 Smm 以 下 。 各 年 4 月 的 莱 散 值 处 


于 较 低 水 平日 燕 散 量 在 0-4.67 mm d", 日 蒸 散 量 均值 为 1.89 mm-d', 5—7 月 蒸 散 量 整 体 
时 增加 趋势 ， 日 燕 散 量 均值 为 5.45 mm-d', 1) 9~10 月 蒸 散 量 水 平 降 低 ， 日 均 蒸 散 量 为 1.78 
mm-d', 


3.3 蒸 散 量 影响 因素 

蒸 散 量 主要 受 区 域 下 热 面 性 质 和 气象 要 素 的 影响 60， 在 其 驱动 下 呈现 出 独特 的 时 空 变 
化 特征 。 
33.1 土地 利用 / 履 被 类 型 对 蒸 散 量 的 影响 


为 进一步 探究 土地 利用 / 履 被 类 型 对 蒸 散 的 影响 ， 以 2015 年 为 例 ， 统 计 该 年 不 同 土地 利 
用 / 覆 被 类 型 的 日 蒸 散 量 像 元 频率 分 布 ， 分 析 不 同 土地 利用 / 履 被 面积 变化 所 导致 的 日 总 蒸 散 
量 的 变化 情况 (图 7)。 水 体 及 河谷 湿地 的 日 蒸 散曲 线 比 较 靠 右 ， 蒸 散 量 较 大 ， 频 率 分 布 随 
着 不 同 的 季节 有 明显 的 变化 。 其 次 为 耕地 ,夏季 耕地 的 日 蒸 散 量 与 河谷 湿地 的 蒸 散 量 较为 接 
近 , 在 2-6 mm 波动 。 草 地 的 蒸 散 量 在 0.5~4 mm 变化 , 不 同 季节 草地 的 频率 分 布 大 致 相同 。 
沙 地 和 农村 居民 点 的 日 蒸 散 量 最 小 , 且 像 元 频率 分 布 比较 集中 , 植被 稀 朴 , 涵养 水 源 能 力 差 ， 
使 得 蒸 散 量 活动 受 限 ， 日 蒸 散 量 在 0-2 mm。 总 体 来 看 ， 植 被 区 蒸 散 量 大 于 非 植被 区 的 蒸 散 
量 ， 蒸 散 主要 与 植被 覆盖 度 和 水 分 的 供给 条 件 有 关 呈 1。 蒸 散 量 依次 为 : 水 体 > 沼泽 地 > 耕地 > 
草地 > 沙 地 > 农村 居民 点 。 


I 


(a)4 月 7 


ET/mm 
一 一 耕地 一 一 农村 居民 点 一 一 沙 地 一 一 水 体 沼泽 地 一 一 草地 
7 2015 年 不 同 土地 利用 类 型 ET 值 统计 对 比 
Fig. 7 Statistical comparison of ET values of different land use types in 2015 

不 同 土地 利用 / 履 被 类 型 导致 ET 差异 的 原因 有 : 植被 覆盖 度 及 植被 类 型 不 同 ，4 月 植被 
返青 ，9 月 植被 到 了 生长 后 期 开始 枯萎 ， 加 之 流域 大 部 分 地 区 牧民 开始 打 草 ， 植 被 覆盖 度 降 
低 ， 流 域 整 体 蒸 散 量 较 小 。 生 长 季 沼 译 植物 密集 ， 植 被 覆盖 度 高 ， 植 被 的 蒸腾 作用 大 ， 沼 泽 
显 地 蒸 散 量 就 越 大 , 这 与 沈 欣 等 中 在 白洋淀 湿地 生态 系统 的 研究 结果 一 致 ; 土壤 类 型 不 同 ， 
究 区 草地 土壤 类 型 主要 有 黑 钙 土 、 厚 栗 黄 土 等 , 沙 地 主要 为 襄 漠 风沙 土 , 沼泽 地 主要 为 革 
原 沼泽 土 ， 不 同 的 土壤 类 型 对 于 太阳 辐射 的 吸收 不 同 中 ， 土 壤 含 水 量 不同 ， 沼 泽 地 有 锡 林 


oe 


k= 


i= mal 
HFI 
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河 补给 ,耕地 基本 配备 喷灌 设施 ， 水 分 供给 充足 , 在 相同 的 条 件 下 沼泽 地 和 耕地 的 蒸发 量 要 
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大 于 研究 区 其 他 土地 利用 类 型 六 。 农 村 居民 点 和 沙 地 的 燕 散 量 较 低 ， 一 般 农村 居民 点 因 其 
包含 大 量 不 透水 层 ， 导 致 地 表 水 源 补给 减少 ， 日 蒸 散 作用 也 相应 减弱 ， 表 现 为 " 干 岛 "效应 ， 


其 通过 吸收 太阳 辐射 能 量 使 地 表 温 度 增加 ， 造 成 感 热 高 潜 热 低 ， 因 此 蒸 散发 量 较 低 占 。 沙 


地 相 比 于 农村 居民 点 而 言 ，ET 一 部 分 来 自 于 表层 土壤 水 及 降水 蒸发 ， 还 有 一 部 分 来 自 于 低 


覆盖 的 沙 地 灌木 从 、 


披 碱 草 等 履 沙 植被 的 蒸腾 ， 所 以 ，ET7 值 较 农 村 居民 点 略 高 。 


3.32 气象 因子 对 蒸 散 量 的 影响 
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8 日 蒸 散 量 与 气象 因子 的 关系 


Fig. 8 The relationship between daily evapotranspiration and meteorological factors 


图 8 是 ET 与 平均 温度 T、VPD、R。、RH 之 间 的 关系 。 可 以 看 出 ，ET 与 Rs 的 拟 合 效果 


最 好 ， 决 定 系数 为 0.83， 其 次 为 T 和 VPD， 决 定 系数 分 别 为 0.71 和 0.63， 而 对 RH 的 响应 


程度 最 低 ， 为 0.26 


o 


由 表 3 AYA T. VPD 和 R E ET 呈正 相关 关系 ， 分 别 为 0.76、0.72、0.83。 与 RH ffi 
相关 关系 ， 为 0.53。 从 敏感 系数 绝对 值 大 小 可 看 出 ，ET 对 ,的 变化 最 敏感 ，REH KZ, T 
和 VPD 的 敏感 性 较 低 。 综 合 来 看 Ry 是 影响 区 域 蒸 散 量变 化 的 的 主要 因子 。 


表 3 气象 要 素 敏感 系数 及 其 与 ET 的 相关 系数 


Table2 Meteorological element sensitivity coefficient and correlation coefficient with ET 


气象 要 素 与 ET 相关 系数 敏感 系数 
平均 温度 (7) 0.76” 0.31 
饱和 水 汽 压 差 (VPD) 0.72™ 0.29 
地 表 净 辐射 (Ru) 0.837 0.76 
相对 湿度 (RH) -0.53™ -0.47 


TE: ** 表 示 在 0.01 性 水 平 ( 双 侧 ) 上 显著 相关 。 


从 根本 上 来 看 


,无 论 是 植物 蒸腾 还 是 土壤 蒸发 , 在 植被 体内 水 分 和 土壤 水 分 发 生 相 变 汽 


化 的 过 程 中 ,太阳 辐 射 是 主要 的 能 量 来 源 ， 太 阳 辐 射 增加 ， 导 致 气温 和 水 温 升 高 ， 水 分 蒸发 


速率 增加 。 因 此 ， 


日 燕 散 量 与 太阳 辐射 存在 较 好 的 线性 关系 ， 这 与 牛 忠 恩 等 "研究 中 国 陆 


uc 


生态 系统 蒸 散 时 空 变 化 及 其 影响 因素 所 得 到 的 结论 一 致 。 沼 译 湿 地 蒸 散 量 与 饱和 水 汽 压 呈 


=i 


lK 


it 


E 相 关 关 系 ， 但 在 饱和 水 汽 压 过 大 之 后 ， 饱 和 水 汽 压 的 增加 反而 抑制 了 蒸 散 量 的 增 大 ， 因 为 
$ 饱 和 水 汽 压 很 大 时 ， 蒸 腾 作 用 强烈 ， 引 起 气孔 部 分 关闭 ， 降 低 燕 腾 失 水 ，Liu 等 中 在 分 析 
E 京 地 区 气候 变化 对 燕 散 量 的 影响 得 到 相同 的 结论 。 空 气相 对 湿度 大 ,水 汽 落 发 强度 小 ， 这 


与 书 振 锋 等 中 研究 陕西 省 蒸 散 量 对 气候 因子 的 响应 所 得 到 的 结果 相同 。 区 域 蒸 散 量变 化 与 
驱动 因子 的 关系 十 分 复杂 , 定量 化 度量 相关 因子 对 流域 蒸 散 的 影响 会 在 进一步 的 研究 中 加 以 


算 锡 林 河 流域 上 游 的 燕 散 量具 有 可 行 性 。 


完善 
4 结论 


o 


CD 利用 SEBAL 模型 估算 的 尺 、G 及 本 值 与 实测 值 较 接 近 ， 故 采用 SEBAL 模型 估 


(2) Rs 和 ET 具有 较 高 的 一 致 性 ， 季 节 性 变化 显著 ， 空 间 上 呈现 东 高 西 低 的 趋势 ， 河 


谷 湿地 的 R 和 ET 在 生长 季 内 保持 较 高 的 水 平 ，R, 值 在 400-750 Wm’ 之 间 ，ET 值 在 
2.42~10.67 mm:d 1， 沙 丘 裸 地 燕 散 量 较 少 ， 在 5 mm UF. 


(3) 土地 利用 / 履 被 类 型 和 气象 因子 是 影响 蒸 散 量 时 空 变化 的 主要 因素 。 不 同 土地 利用 


类 型 蒸 散 量 依次 为 : 水 体 > 沼 泽 地 > 耕地 > 草地 > 沙 地 > 农村 居民 点 ; 锡 林 河 流域 上 游 ET 5 T. 
VPD, R, 呈正 相关 ， 与 RH 时 负 相 关 。 从 敏感 系数 来 看 ，ET 对 Ri 的 变化 最 敏感 。 
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Estimation of evapotranspiration and its spatiotemporal 
characteristics in the upper reaches of the 
Xilin River Basin 
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Abstract: To reveal the evapotranspiration and its temporal and spatial distribution characteristics in the upper 
reaches of the Xilin River Basin, The ET estimation of 12 growth seasons during 2015—2017 was conducted using 
a Surface Energy Balance Algorithm for Land (SEBAL) model with Landsat 8 data and meteorological 
observations, and analysis of its temporal and spatial patterns and its influencing factors.The results show that :(1) 
The inversion parameter values of the model are well fitted to the measured values, the decision coefficients of 
surface net radiation , soil heat flux, ET are 0.80, 0.65, 0.86.(2) The spatial distribution of R, and ET at the same 
period has a high consistency, and the overall trend is east-high and low-west, with significant spatial 
differentiation characteristics. The daily ET value in April is 0~4.67 mm, from May to July, the daily ET increased 
overall, the maximum value is 10.32mm, it's decreased from 0 to 3.34 mm from September to October. The 
evapotranspiration of different land-use types is: water > swamp > arable land > grassland > sandy land > rural 
settlements.(3) The evapotranspiration is positively correlated with the average temperature, saturated vapor 
pressure, and solar net radiation, negatively correlated with relative humidity, ET is most sensitive to changes in 


solar radiation. 


Key words: SEBAL model; Evaporation; The space-time characteristics; Xilin River Basin 


